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In SU(3) lattice QCD, we improve a method to extract gauge fields from link variables analyti- 
cally. With this method, we perform the first study on the effective mass generation of off-diagonal 
gluons and infrared Abelian dominance in the maximally Abelian (MA) gauge in the SU(3) case. 
We investigate the propagator and the effective mass of the off-diagonal gluon field in the MA 
gauge with U(l) 3 x U(l) g Landau gauge fixing in SU(3) quenched lattice QCD on 16 4 at j3=5.7, 
5.8 and 6.0. The off-diagonal gluon component behaves as a massive vector boson with the 
approximate effective mass M s = 1.1 — 1 .2GeV in the region of r = 0.3 — 0.8fm, and its propa- 
gation is limited within a short range. We thus show the origin of infrared Abelian dominance in 
terms of short-range propagation of off-diagonal gluons. We also investigate the functional form 
of the off-diagonal gluon propagator. We find that the functional form is well described by the 
four-dimensional Euclidean Yukawa-type function exp{— m ffr} / r with m ff = 1 .3 — 1 .4GeV for 
r = 0.1 — 0.8fm. This also indicates that the spectral function has a negative region. 
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1. Introduction 



For the quark-confinement mechanism, Nambu, 't Hooft and Mandelstam proposed the dual- 
superconductor picture [JJ. This picture is based on the electromagnetic duality and the analogy 
with the one-dimensional squeezing of the magnetic flux in the type-II superconductor. In this pic- 
ture, there occurs color-magnetic monopole condensation, and then the color-electric flux between 
the quark and the antiquark is squeezed as a one-dimensional tube due to the dual Higgs mech- 
anism. From the viewpoint of the dual-superconductor picture in QCD, however, there are two 
assumptions of Abelian dominance [Q, |5|] and monopole condensation. Here, Abelian dominance 
means that only the diagonal gluon component plays the dominant role for the nonperturbative 
QCD phenomena like confinement. 

The maximally abelian (MA) gauge has mainly been investigated from the viewpoint of the 
dual-superconductor picture [§ 1 1| 0, | £§ ^Oj, [IT} [T| 0, [H] and the various lattice QCD 
Monte Carlo simulations show that the MA gauge fixing seems to support these assumptions [Q, D, 

According to these studies, the diagonal gluons seem to be significant to the infrared QCD 
physics, which is called "infrared Abelian dominance". Infrared Abelian dominance means that 
off-diagonal gluons do not contribute to infrared QCD. Therefore, the essence of infrared Abelian 
dominance is the behavior of the off-diagonal gluon propagator. 

The gluon propagators in the MA gauge has been investigated in SU(2) lattice Monte Carlo 



simulations [11, 14, 16]. To investigate the gluon propagators in the MA gauge, it is desired to ex- 
tract the gluons exactly from the link- variables, because the link- variable cannot be expanded even 
for a small lattice spacing due to large fluctuation of gluons. In SU(2) lattice case, the extraction 
is easy to be done without any approximation, because of the SU(2) property. With this extraction, 
the SU(2) lattice simulation suggests that the off-diagonal gluons do not propagate in the infrared 
region due to the effective mass M g- ~ 1.2GeV, while the diagonal gluon widely propagates [11]. 



In this paper, we propose a method to extract the gluons from the link-variable directly and 
generally in SU(3) lattice QCD, and to investigate the gluon propagators in the MA gauge. 



2. SU(3) lattice QCD results for gluon propagators in the MA gauge 



To begin with, we consider a useful and general method to extract the gauge fields analytically 



and exactly from the link- variables whether |agA^(3c)| 1 is satisfied or not [17, 18]. 

In this method, each link- variable U = (s) is diagonalized with a unitary matrix £2, 



U d = ClUtf = e iagQAa1 = 



/e' e > \ 



,102 



V 







J&3 



(2.1) 



where — n < 0, < K {i = 1,2,3) is taken. Note that 6\ + 02 + 6s = (mod 2ti) according to 
TrfiAQ' = TrA = (mod 2k). This property is also numerically checked in the MA gauge with 
U(l)3xU(l)8 Landau gauge fixing. 
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We can derive gluon fields A by taking the logarithm of U ( /, 



QA& = — 

ag 



( 0i \ 




(2.2) 



This formalism is quite general, because the derivation is correct in any gauge and even without 
any gauge fixing. 



Using the SU(3) lattice QCD, we calculate the gluon propagators Q18Q in the MA gauge with 
the U(l)3xU(l)8 Landau gauge fixing. In the MA gauge, to investigate the gluon propagators, 
we use the above-mentioned method. The Monte Carlo simulation is performed with the standard 
plaquette action on the 16 4 lattice with /3 =5.7, 5.8 and 6.0 at the quenched level. All measurements 
are done every 500 sweeps after a thermalization of 10,000 sweeps using the pseudo heat-bath 
algorithm. We prepare 50 gauge configurations for the calculation at each /3. The statistical error 
is estimated with the jackknife method. 

Here, we study the Euclidean scalar combination of the diagonal (Abelian) and off-diagonal 
gluon propagators as 

L o=3,8 

The scalar combination of the propagator is expressed as the function of the four-dimensional 
Euclidean distance r = (xn — y^) 2 - When we consider the renormalization, these propagators 
are multiplied by an r-independent constant, according to a constant renormalization factor of the 
renormalized gluon fields. 

We show in Fig.[l| the lattice QCD result for the diagonal gluon propagator G^jf (r) and the off- 
diagonal gluon propagator G^(r) in the MA gauge with the U(l)3xU(l)8 Landau gauge fixing. 
In the MA gauge, G^ l (r) and G?S(r) manifestly differ. The diagonal-gluon propagator G^ l (r) 
takes a large value even at the long distance. In fact, the diagonal gluons in the MA gauge 

propagate over the long distance. In contrast, the off-diagonal gluon propagator (r) rapidly 
decreases and is negligible for r > 0.4fm in comparison with G^ l (r). Then, the off-diagonal 
gluons A" (a ^ 3,8) seem to propagate only within the short range as r < 0.4fm. Thus, "infrared 
abelian dominance" is found in the MA gauge. 

3. Estimation of off-diagonal gluon mass in the MA gauge 



Next, we investigate the effective mass of off-diagonal gluons [18]. We start from the La- 
grangian of the free massive vector field A^ with the mass M / in the Euclidean metric. In the 
infrared region with large Mr, the propagator G^^(r;M) reduces to 

3-y/M e~ Mr 

G m {r,M) = (A^A^y)) ~ / - 3 , (3.1) 

2(271) 2 ri 
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Figure 1: The SU(3) lattice QCD results of the gluon propagators G^Jf'fr) and G°^(r) as the function of 
r = J (xfi —y/i) 2 in the MA gauge with the U(l)3xU(l)g Landau gauge fixing in the physical unit. The 
Monte Carlo simulation is performed on the 16 4 lattice with /3 = 5.7, 5.8 and 6.0. The diagonal-gluon 
propagator G^ el (r) takes a large value even at the long distance. On the other hand, the off-diagonal gluon 
propagator G?E (r) rapidly decreases. 




r [fm] 

Figure 2: The logarithmic plot of r^'^G^ (r) and r 3 / 2 G^ el (r) as the function of the Euclidean distance r 
in the MA gauge with the U(l)3xU(l)g Landau gauge fixing, in the SU(3) lattice QCD with 16 4 at j3 = 5.7, 
5.8 and 6.0. The solid line denotes the logarithmic plot of r 3 / 2 G^ (r) ~ r l l 2 K\ {Mr) in the Proca formalism. 



In Fig^, we show the logarithmic plot of r 3 ' 2 G^(r) and r^l 2 G^jf(r) as the function of the 
Euclidean distance r in the MA gauge with the U(l)3 xU(l)g Landau gauge fixing. From the linear 
slope on r^' 2 G au{r) in the range of r = 0.3 — 0.8 fm, the effective off-diagonal gluon mass M ff is 
estimated. Note that the gluon-field renormalization does not affect the gluon mass estimate, since 
it gives only an overall constant factor for the propagator. We summarize in Table 1 the effective 
off-diagonal gluon mass M g obtained from the slope analysis at j8 =5.7, 5.8 and 6.0. Therefore, 
the off-diagonal gluons seem to have a large mass M g- ~ 1 . 1 — 1 .2 GeV. This result approximately 
coincides with SU(2) lattice calculation [11]. 

Finally in this section, we discuss the relation between infrared abelian dominance and the 
off-diagonal gluon mass. Due to the large effective mass M ff, the off-diagonal gluon propagation 
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Table 1: Summary table of conditions and results in SU(3) lattice QCD. In the MA gauge, the off-diagonal 
gluons seem to have a large effective mass M ff ~ 1.1 — 1.2 GeV and the functional form in the range of 
r = 0.1 — 0.8 fm is well described with the four-dimensional Euclidean Yukawa function ~ exp(— m ^r)/r 
with OToff ~ 1.3 — 1.4 GeV. 



lattice size 


P 


a[fm] 


M off [GeV] 


m off [GeV] 




5.7 


0.186 


1.2 


1.3 


16 4 


5.8 


0.152 


1.1 


1.3 




6.0 


0.104 


1.1 


1.4 



is restricted within about M~ s ~ 0.2fm in the MA gauge. Therefore, at the infrared scale as r> 
0.2fm, the off-diagonal gluons A^ (a ^ 3, 8) cannot mediate the long-range force like the massive 
weak bosons in the Weinberg-Salam model, and only the diagonal gluons A 3 , A^ can mediate the 
long-range interaction in the MA gauge. In fact, in the MA gauge, the off-diagonal gluons are 
expected to be inactive due to the large mass M e in the infrared region in comparison with the 
diagonal gluons. Then, infrared abelian dominance holds for r 3> . 

4. Analysis of the functional form of off-diagonal gluon propagator in the MA gauge 

In this section, we investigate the functional form of the off-diagonal gluon propagator G°Jl (r) 



in the MA gauge in SU(3) lattice QCD [18]. In the previous section, we compare the off-diagonal 
gluon propagator with the massive vector boson propagator and estimate the gluon mass. In fact, 
the gluon propagator would not be described by a simple massive propagator in the whole region 
of r = 0.1 -0.8 fin. 



There is the similar situation in the Landau gauge Q19[|. The functional form of the gluon 
propagator cannot be described by exp(— Mr)/r 3 / 2 with an effective mass M in the whole region 
of r = 0.1 — 1.0 fm. The appropriate form is the four-dimensional Euclidean Yukawa-type function 
exp(— mr)/r with a mass parameter m. 

In the same way, in the MA gauge, we also compare the gluon propagator with the four- 
dimensional Euclidean Yukawa function. In Fig.|3[ we show the logarithmic plot of rG^(r) and 

^Abel 



rGuu'fr) as the function of the distance r in the MA gauge with the U(l)3xU(l)8 Landau gauge 



fixing. Note that the logarithmic plot of rG^Ar) is almost linear in the whole region of r = 0.1 — 
0.8 fm, and therefore the off-diagonal gluon propagator is well expressed by the four-dimensional 
Euclidean Yukawa function Ae~ m ° Br /r, with a mass parameter m D ff and a dimensionless constant 
A. The best-fit mass parameter m g- is given in Table 1 at each j8 = 5.7, 5.8 and 6.0. 



We comment on the four-dimensional Euclidean Yukawa-type propagator Q19|]. If the func 



tional form of the off-diagonal gluon is well described by the four-dimensional Yukawa function, 
we analytically calculate the off-diagonal zero-spatial-momentum propagator, D^^) = J d 3 xG° t ^ l (r), 
and obtain the spectral function by the inverse Laplace transformation. Similarly in the Landau 
gauge [|TJ], we thus derive the spectral function p oS (co) of off-diagonal gluons in the MA gauge, 

p ° tt ((o) = -— rw 0(w - m off + — -£ —8 (co - woff). (4.1) 
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Figure 3: The logarithmic plot of rG^(r) and rG^j~\r) as the function of the Euclidean distance r in the 
MA gauge with the U(l)3xU(l)g Landau gauge fixing, using the SU(3) lattice QCD with 16 4 at J3=5.7, 5.8 
and 6.0. For rG?E (r), the approximate linear correlaton is found. 



5. Summary and Concluding Remarks 

We have performed the first study of the gluon propagators in the MA gauge with the U(l)3 xU(l)g 
Landau gauge fixing in the SU(3) quenched lattice QCD. To investigate the gluon propagators in 
the MA gauge, we have considered to derive the gluon fields analytically from the SU(3) link- 
variables. 

With this method, we have calculated the Euclidean scalar combination G^^(r) of the diagonal 
and the off-diagonal gluon propagators, and have considered the origin of infrared Abelian dom- 
inance. The Monte Carlo simulation is performed on the 16 4 lattice at /3=5.7, 5.8 and 6.0 at the 
quenched level. We have found that the off-diagonal gluons behave as massive vector bosons with 
the effective mass M g ~ 1.1 — 1.2 GeV for r = 0.3 — 0.8 fm. The effective gluon mass has been 
estimated from the linear fit analysis of the logarithmic plot of r 3 / 2 G?S(r). Due to the large value, 
the finite-size effect for the off-diagonal gluon mass is expected to be ignored. The large gluon 
mass shows that the off-diagonal gluons cannot mediate the interaction over the large distance as 
r 3> M~ff , and such an infrared inactivity of the off-diagonal gluons would lead infrared Abelian 
dominance in the MA gauge. 

On the other hand, from the behavior of the diagonal gluon propagator G^ el (r) and r 3 / 2 G^ el (r), 



the diagonal gluons seem to behave as light vector bosons []18j]. For the detailed argument on 



G^jf'(r), one should consider the finite size effect more carefully, because the diagonal gluons 
would propagate over the long distance beyond the lattice size. 

Finally, we have also investigated the functional form of the off-diagonal gluon propagator 
G^(r) in the MA gauge. We show that G^(r) is well described by the four-dimensional Eu- 
clidean Yukawa-type form with the mass parameter m fr ^ 1.3 — 1.4 GeV in the whole region of 
r = 0.1 — 0.8 fm. This indicates that the spectral function p off (ft)) of the off-diagonal gluons in the 
MA gauge has the negative- value region [18], as in the Landau gauge [|l^, 20, 21]. 
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In this study, we investigate the off-diagonal gluon propagator. To be strict, the off-diagonal 
gluon propagator consists of two scalar functions corresponding to longitudinal and transverse 
components. Therefore, we will investigate each effective mass and the functional form of these 
components. 
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